Robustness of quantum discord to sudden death 
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We calculate the dissipative dynamics of two-qubit quantum discord under Markovian environ- 
ments. We analyze various dissipative channels such as dephasing, depolarizing, and generalized 
amplitude damping, assuming independent perturbation, in which each qubit is coupled to its own 
channel. Choosing initial conditions that manifest the so-called sudden death of entanglement, we 
compare the dynamics of entanglement with that of quantum discord. We show that in all cases 
where entanglement suddenly disappears, quantum discord vanishes only in the asymptotic limit, 
behaving similarly to individual decoherence of the qubits, even at finite temperatures. Hence, 
quantum discord is more robust than the entanglement against to decoherence so that quantum 
algorithms based only on quantum discord correlations may be more robust than those based on 
entanglement. 
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Entanglement is widely seen as the main reason for the 
computational advantage of quantum over classical algo- 
rithms. This view is backed up by the discovery that, 
in order to offer any speedup over a classical computer, 
the universal pure-state quantum computer would have 
to generate a large amount of entanglement [T. However, 
quantum entanglement is not necessary for deterministic 
quantum computation with one pure qubit (DQCl), in- 
troduced by Knill and Laflamme 2 . As seen in [3 , 4 , al- 
though there is no entanglement, other kinds of nonclas- 
sical correlation are responsible for the quantum compu- 
tational efficiency of DQCl. Such correlations are char- 
acterized as quantum discord which accounts for all 
nonclassical correlations present in a system, being the 
entanglement a particular case of it. Besides its appli- 
cation in DQCl, quantum discord has also been used 
in studies of quantum phase transition 0, estimation 
of quantum correlations in the Grover search algorithm 
[7 and to define the class of initial system-bath states 
for which the quantum dynamics is equivalent to a com- 
pletely positive map [8J. 

When considering a pair of entangled qubits exposed 
to local noisy environments, disentanglement can occur 
in a finite time UHl HIl [H US], differently from the 
usual local decoherence in asymptotic time. The occur- 
rence of this phenomenon, named "entanglement sudden 
death" (ESD), depends on the system-environment in- 
teraction and on the initial state of the two qubits. Our 
goal was to investigate the dynamics of quantum discord 
of two qubits under the same conditions in which ESD 
can occur. We show in this letter that even in cases where 
entanglement suddenly disappears, quantum discord de- 
cays only in asymptotic time. Furthermore, this occurs 
even at finite temperatures. In this sense, quantum dis- 
cord is more robust against decoherence than entangle- 



ment, implying that quantum algorithms based only on 
quantum discord correlations are more robust than those 
based on entanglement. 

There are various methods to quantify the entangle- 
ment between two qubits [Ml [T5j [1^ and, even when 
they give different results for the degree of entanglement 
of a specific state, all of them result in for separable 
states. Therefore, under dissipative dynamics where an 
initial entangled state can disappear suddenly, all of them 
necessarily agree about the time when the quantum state 
becomes separable. 

Here, to investigate the two-qubit entanglement dy- 
namics we use concurrence as the quantifier [15 . The 
concurrence is given by max {0, A (t)}, where A(t) = Ai — 
A2 — A3 — A4 and Ai > A2 > A3 > A4 are the square roots 
of the eigenvalues of the matrix p{t)a2 <8) <T2P*{t)a2 'S) 0^2, 
p*{t) being the complex conjugate of p{t) and a2 the sec- 
ond Pauli matrix. The density matrix we use to eval- 
uate concurrence has an X structure [12j . defined by 
P12 — P13 — P24 = P34 — 0, which are constant during the 
evolution, for the various dissipative channels used here. 
In this case, the concurrence has a simple analytic ex- 
pression C (t) = 2 max {0, Ai (t) , A2 (t)} , where Ai (t) = 
IpwI - yJpi2Pzz and A2 (i) = |p23| - V^iipil- However, 
as pointed out above, entanglement is not the only kind 
of quantum correlation. In quantum information theory, 
the Von Neumann entropy, S'(p) — — Tr (plog p), is used 
to quantify the information in a generic quantum state p. 
The total correlation between two subsystems A and B of 
a bipartite quantum system is given by the mutual 
information, 

'^{pab)^S{pj^ + S{pb)-S{pab)-. (1) 

where ^(/O^g) = — Tr (p^g log /O^g) is the joint entropy 
of the system [T71. This bipartite quantum state is 
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a hybrid object with both classical and quantum char- 
acteristics and, in order to reveal the classical aspect of 
correlation, Henderson and Vedral suggested that corre- 
lation could also be split into two parts, the quantum and 
the classical [TH]. The classical part was defined as the 
maximum information about one subsystem that can be 
obtained by performing a measurement on the other sub- 
system. If we choose a complete set of projectors {11^} 
to measure one of the subsystems, say B, the information 
obtained about A after the measurement resulting in out- 
come k with probability pk , is the difference between the 
initial and the conditional entropy [T^ 

Qa{pab) = max[S'(p^) - Vpfc5'(p^|fc)], (2) 
{n,} ^ 

where = Tr^iUkPAB'^k) /T^rABO^kPABT^k) is the 

reduced state of A after obtaining the outcome k 
in B. This measurement of classical correlation as- 
sumes equal values, irrespective of whether the mea- 
surement is performed on the subsystem A or B, 
for all states p^B such that S{pa) ~ S{pb) 
This condition is true of all density operators used 
in this paper since they can be written as p — 

\ I + Co ((73 (g) I + I (Ts) + Ylij ^■j'^j ® ^3 1 where Ci and 
ai {i — 1,2,3) are real constants and Pauli matrices, re- 
spectively. Therefore pa = PB- 

In this scenario, a quantity that provides information 
on the quantum component of the correlation between 
two systems can be introduced as the difference between 
the total correlations in ([l| and the classical correlation 
in (j2]). This quantity is identical to the definition of 
quantum discord introduced by OUivier and Zurek in [5J , 
namely 



D (pab) = I (pab) - Q (pab) , 



(3) 



this being zero for states with only classical correlations 
[51 [TH] and nonzero for states with quantum correlations. 
Moreover, quantum discord includes quantum correla- 
tions that can be present in states that are not entan- 
gled [S], revealing that all the entanglement measure- 
ments such as concurrence, entanglement of formation, 
etc, do not capture the whole of quantum correlation be- 
tween two mixed separate systems. For pure states, the 
discord reduces exactly to a measure of entanglement, 
namely the entropy of entanglement. 

In order to calculate the quantum discord between two 
qubits subject to dissipative processes, we use the follow- 
ing approach. The dynamics of two qubits interacting in- 
dependently with individual environments is described by 
the solutions of the appropriate Born-Markov-Lindblad 
equations [19j , that can be obtained conveniently by the 
Kraus operator approach \17\ t ] . Given an initial state for 
two qubits p(0), its evolution can be written compactly 
as 



where the so-called Kraus operators = E^^E^ [T7] 
satisfy T,^^^Ej^ j,Efj,^i, = I for all t. The operators 
describe the one-qubit quantum channel effects. 

In the cases where the quantum channel induces a dis- 
entanglement only in asymptotic time, the quantum dis- 
cord does not disappear in a finite time, since the entan- 
glement is itself a kind of quantum correlation. There- 
fore, we present below what happens to the discord in 
the ESD situations, for some of the common channels for 
qubits: dephasing, generalized amplitude damping (ther- 
mal bath at arbitrary temperature) and depolarizing. 

Dephasing: The dephasing channel induces a loss 
of quantum coherence without any energy exchange 
|17j . The quantum state populations remain unchanged 
throughout the time. To examine the two-qubit entan- 
glement and discord dynamics under the action of a de- 
phasing channel, we utilize the Werner state as the initial 
condition, i.e, p{0) = (1 - a)I/4 + a l^*") a G [0, 1] 

and \^^) = (|01) - |10)) /\/2. In this case, we can cal- 
culate analytically the quantum discord in a situation 
where the entanglement suddenly disappears. Thus, ac- 
cording to Eq. Q, with the non-zero Kraus operators 
for a dephasing channel given by Eq = diag (l, \/r^^) 
and El ~ diag (l, -^7), where 7=1 — e^'"*, T denoting 
the decay rate |17j . the elements of the density matrix of 
this system evolve to 



Pti{t) 

P23{t) 



pu{0),i= 1..4, 
P23(0)(l-7) = 



P32(i)- 



The concurrence for this state is given by C(p) = 
a (3/2 — 27) — 1/2, which reaches zero in a finite time 
for any a 7^ 1, as shown in Fig. 1(a). On the other hand, 
based on the results given in [20!, the quantum discord 
for this state reads D (p) = [F{a + b) + F{a - b)] /4 - 
F{a)/2, where F{x) — xlogjX, a = (1 — a) and b = 
2a(l — 7). As shown in Fig. 1(b), for any a, the quan- 
tum discord vanishes {D (p) = 0) only in the asymptotic 
limit. 




Pit) 



^^,,.E^,,p{0)El,, 



(4) 



FIG. 1: Dissipative dynamics of (a) concurrence and (b) dis- 
cord as functions of a and 7, assuming independent dephasing 
perturbative channels. 

Generalized Amplitude Damping (GAD): The 
GAD describes the exchange of energy between 
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the system and the environment, including fi- 
nite temperature aspects. It is described by the 
Kraus operators Eq — y^diag (1,^1 ~ 7)' — 
^{ai+iG2)/2, E2 = v/(l- 9) 7diag(yr^,l), 
and — ■\/(l — (7) 7 {<7i — i<T2) /2, where 7 is defined 
above and q defines the final probability distribution of 
the qubit when t ^ 00 {q — \ corresponds to the usual 
amplitude damping with T = QK) |17j . 

For the initial condition given by p (0) = |<i>) (<i>| with 

1$) = Vr^lOO) + Va|ll) , ae[0,l], (5) 

we obtain, according to Eq.Q, the density matrix dy- 
namics 

Pii(i) = pii(0){l-7[2(l--7)-7(l-2g)]} + 7V, 
P22 (t) = P33 (t) = 7 [pii (0) (1 - 2q) (1 - 7) + g (1 - jq)] , 

Piiit) = 1 - pii (t) - 2p22 (i) , 

PU (t) = PAl (t) = P14 (0) (1 - 7) • 

We examine the dissipative dynamics derived from this 
channel, taking q — 1 and q — 2/3. For these cases, we 
compute the discord numerically and compare it with 
the concurrence. To calculate the discord, we chose 
the set of projectors {|V'i) {i'll , 1^2) {ip2\}, where IV'i) = 
cos \0)+e''l' sin 9 \ 1) and | V'2) = - cos \ 1) -t-e""^ sin 0\O), 
to measure one of the subsystems. The maximum of 
equation ([2]) is obtained numerically by varying the an- 
gles and 4> from to 2tt. 

For q = 1, we have C{t) = max{0, Ai(i)}, with 
Ai{t) > for all t, whenever y/1 — a > y^,see Fig. 2(a). 
On the other hand, for q ^ 1, we have Ai(t — ^ 00) < 
for all a, as shown in Fig. 2(c). In this case, since the 
concurrence decays monotonically under the Markovian 
approximation [5T] , the ESD occurs for any initial super- 
position state, i.e, for any a different from or 1. 

The discord evidently behaves differently from the con- 
currence, see Fig. 2(b) and 2(d). In both situations 
(g = 1 and q = 2/3), the discord decays exponentially 
and vanishes only asymptotically. For a = 1 or a = 
(pure separable states) the quantum discord is zero all 
the time, as expected. 

Depolarizing: The depolarizing channel represents the 
process in which the density matrix is dynamically re- 
placed by the maximal mixed state 1/2, I being the iden- 
tity matrix of a single qubit. The set of Kraus op- 
erators that reproduces the effect of the depolarizing 
channel is given by Eq = ^1 — 37/41, Ei = ^/jjiax, 
E2 = \/7/4<7y, and E^ — y/^/Aaz, with 7 as defined 
above |T7]. Assuming the initial condition given by Eq. 
([5]), we obtain the density matrix element dynamics: 

Pii(i) - pii(0)(l-7)+7V4, 

P22{t) = P33(i) =7/2(1-7/2), 

P4i{t) = l-pii(t)-2p22(i), 

pu{t) = p4i(i) =Pi4(0)(l-7), 




FIG. 2; Dissipative dynamics of concurrence (a,c) and discord 
(b,d) as a functions of a and 7, assuming independent gener- 
alized amplitude damping. (a,b) q = 1 and (c,d) q = 2/3. 

As in generalized amplitude damping, the ESD occurs 
for any initial condition, since C{t) — max{0, Ai(t)} and 
Ai{t — > 00) < 0, Fig. 3(a). Again, as shown in Fig. 3(b), 
the quantum discord does not disappear in a finite time. 
Here we used the same procedure as above to calculate 
the discord numerically. 




FIG. 3: Dissipative dynamics of (a) concurrence and (b) dis- 
cord as functions of a and 7, assuming independent depolar- 
izing channels. 

Dephasing plus Amplitude Damping: In Fig 4 we 
plot the concurrence and quantum discord for the case 
where both qubits interact individually with two distinct 
reservoirs: those which induce dephasing and amplitude 
damping. We assume equal decay rates (F) for both 
channels, q — 1 {T — OK) and 7 = 1 — e^'"*. For 
the initial condition given by Eq. ([5|, the dephasing 
channel alone is not able to induce sudden death of en- 
tanglement. On the other hand, in the presence of an 
amplitude damping channel, the entanglement suddenly 
disapears for some values of a, Fig. 2(a), as discussed 
above. However, when both channels are present, the 
dynamics of the entanglement is very different, suddenly 
disappearing for all values of a as we can see in Fig 4(a). 
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This non-additivity of the decoherence channels in the 
entanglement dynamics was firstly pointed out by T. Yu 
and J. H. Eberly [11], in contrast to the additivity of 
de decay rates of different decoherence channels of a sin- 
gle system [TT]. But, as shown in Fig 4(b), the discord 
still decays asymptotically when both decoherence chan- 
nels are present, as shown in Fig 4(b), indicating that 
the additivity of the decoherence channels is valid for the 
quantum discord. 




FIG. 4: Dissipative dynamics of (a) concurrence and (b) dis- 
cord as functions of a and 7, assuming simultaneous action 
of dephasing and generalized amplitude damping channels in 
each qubit with g = 1. 



In conclusion, we have calculated the discord dynamics 
for two qubits coupled to independent Markovian envi- 
ronments. We observed that under the dissipative dy- 
namics considered here, discord is more robust than en- 
tanglement, even at a finite temperature, being immune 
to a "sudden death". This also points to a fact that 
the absence of entanglement does not necessarily indi- 
cate the absence of quantum correlations. Thus, quan- 
tum discord might be a better measure of the quantum 
resources available to quantum information and compu- 
tation processes. This also suggests that quantum com- 
puters based on this kind of quantum correlation, differ- 
ently from those based on entanglement, are more resis- 
tant to external perturbations and, therefore, introduce 
a new hope of implementing an efficient quantum com- 
puter. Moreover, discord may be considered, in this sce- 
nario, as a good indicator of classicaKty [51 [53], since it 
vanishes only in the asymptotic limit, when the coher- 
ence of the individual qubits disappears. However, we 
have not demonstrated that the sudden death of discord 
in a Markovian regime is impossible, and a study of the 
discord from a geometrical point of view |22j , for exam- 
ple, might be useful to address this important question. 
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